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We investigate quantum transport in a normal/superconductor graphene heterostructure, including the possi- 
bility of an anisotropic pairing potential in the superconducting region. We find that under certain circumstances, 
the conductance displays an undamped, oscillatory behaviour as a function of applied bias voltage. Also, we 
investigate how the conductance spectra are affected by a d-wave pairing symmetry. These results combine 
unusual features of the electronic structure of graphene with the unconventional pairing symmetry found for 
instance in high-Tc superconductors. 
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Graphene is a monoatomic layer of graphite with a hon- 
eycomb lattice structure fl']. The electronic properties of 
graphene display several intriguing features, such as a six- 
point Fermi surface and Dirac-like low-energy energy dis- 
persion around the Fermi-points. Condensed matter systems 
with such 'relativistic' electronic structure properties consti- 
tute fascinating examples of low-energy emergent symmetries 
(in this case Lorentz-invariance). Another example where pre- 
cisely this occurs is in one-dimensional interacting fermion 
systems, where phenomena Uke breakdown of Fermi-liquid 
theory and spin-charge separation take place. Graphene fea- 
tures certain similarities to, but also important differences 
from, the nodal Dirac fermions emerging in the low-energy 
sector of the pseudogap phase of d-wave superconductors 
such as the high-Tc cuprates. When Lorentz-invariance 
emerges in the low-energy sector of higher-dimensional con- 
densed matter systems, it is bound to attract much interest 
from a fundamental physics point of view. 

Various aspects of resonant tunneling phenomena in N/N 
and N/I/N graphene structures have recently been investigated 
Although superconductivity does not appear intrinsi- 
cally in graphene, it may nonetheless be induced by means 
of the proximity effect [3]. Motivated by this, the authors of 
Refs. H m considered quantum transport in N/S and N/I/S 
graphene junctions for the case where the pairing potential is 
isotropic, leading to s-wave superconductivity. However, the 
hexagonal symmetry of the graphene lattice also admits un- 
conventional order parameters such as p-wave or d-wave. The 
possible pairing symmetries on a hexagonal lattice up to /- 
wave pairing (I — 3) was given in Ref. |@]. Interestingly, 
among the allowed order parameters, one finds the dj.2_y2- 
symmetry, which is believed to be the dominant pairing sym- 
metry in high-Tc superconductors. Consequently, it should be 
possible to induce superconductivity with nodes in the gap in 
graphene by manufacturing heterostructures of graphene and 
unconventional superconductors. It is of interest to investigate 
how this would affect coherent quantum transport injunctions 
with normal and superconducting graphene. In particular, it 
is essential to study possible zero-energy states (ZES) at the 
interface of such a junction. Such states are known to give rise 



to zero-bias conductance peaks (ZBCPs) in metallic N/S junc- 
tions ff], and will influence the conductance spectra of N/I/S 
junctions. 

In this Letter, we take into account the possibility of an 
anisotropic pairing potential induced in graphene, and study 
coherent quantum transport in both N/S and N/I/S junc- 
tions. In addition, we show that in the latter structure, novel 
conductance-oscillations as a function of bias voltage are 
present both for s-wave and d-wave symmetry of the super- 
conducting condensate due to the presence of low-energy 'rel- 
ativistic' nodal fermions on the N-side. The period of the os- 
cillations decreases with increasing width w of the insulating 
region, and persists even if the Fermi energy in I is strongly 
shifted. This contrasts sharply to metallic N/I/S junctions, 
where the presence of a potential barrier causes the transmit- 
tance of the junction to go to zero with increasing w. The 
feature of conductance-oscillations is thus unique to N/I/S 
junctions with low-energy Dirac-fermion excitations. More- 
over, we contrast the N/S or N/I/S conductance spectra for the 
cases where s-wave and dx2-y2-wa\e superconductor consti- 
tutes the S-side. The former has no nodes in the gap and lacks 
Andreev bound states. The latter has line-nodes that always 
cross the Fermi surface in the gap, and thus features in ad- 
dition to Andreev bound states, also nodal relativistic low- 
energy Dirac fermions. The quantum transport properties in 
a heterostructure of two such widely disparate systems, both 
featuring a particular intriguing emergent low-energy symme- 
try, is of considerable importance. 

The Brillouin zone of graphene is hexagonal and the en- 
ergy bands touch the Fermi level at the edges of this zone, 
amounting to six discrete points. Out of these only two are in- 
equivalent, denoted K and K' and referred to as Dirac points. 
The energy dispersion in the Brillouin zone was calculated 
within a tight-binding model fsl], revealing a conical struc- 
ture of the conduction and valence bands close to the six 
Fermi points, giving rise to an essentially linear dispersion. 
Graphene N/S interfaces contain a new phenomenology com- 
pared to their metallic counterpart, namely the possibility of 
specular Andreev-reflection (AR) fl- In the process of nor- 
mal AR, an incident electron from the N side is reflected as 
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a hole which retraces the trajectory of the electron. In specu- 
lar AR, the reflected hole follows the trajectory which a nor- 
mally reflected electron would have. Depending on whether 
the graphene is doped or not, specular and normal AR will 
compete with each other, also depending on the position of 
the Fermi level with respect to the gap. 

In order to treat the scattering processes at the interfaces 
of the N/I/S junction, we make use of the full Bogoliubov-de 
Gennes (BdG) equation for the 2D sheet of graphene in the 
xy-plane, assuming the clean limit. These equations read 
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where E is the excitation energy, and is the wave-function. 
We use for 4x4 matrices, .?. for 2x2 matrices, and bold- 
face notation for three-dimensional row vectors. Assuming 
that the superconducting region is located at a; > and ne- 
glecting the decay of the order parameter in the vicinity of 
the interface [9], we may write for the spin-singlet order pa- 
rameter Ak = A{9)&'^Q{x), where Q{x) is the Heaviside 
step function, while is the phase corresponding the globally 
broken U{1) symmetry in the superconductor. We consider 
the weak-coupling limit with the momentum k fixed on the 
Fermi surface, such that Ak only has an angular dependence 
9 = atan{ky/kx)- Note that in contrast to previous work, we 
allow for the possibility of unconventional superconductivity 
in the graphene layer since Ak now may be anisotropic. 

Postulating a spin-singlet even parity order parameter, the 
condition A{6) — A{tt + 9) must be fulfilled. The single- 
particle Hamiltonian is given hy H — diag{H^, H^), H± = 
—ivp{axdx±^ydy) . Here, vp is the energy-independent Fermi 
velocity for graphene, while ai denotes the Pauli matrices. 
For later use, we also define the Pauli matrix vector a = 
{ax, <Jy,<Jz)- These Pauli matrices operate on the sublattice 
space of the honeycomb structure, corresponding to the A and 
B atoms, while the ± sign refers to the two so-called valleys 
of K and K' in the Brillouin zone. The spin indices may 
be suppressed since the Hamiltonian is time-reversal invari- 
ant. In addition to the spin degeneracy, there is also a valley 
degeneracy, which effectively allows one to consider either 
the one of the H± set. We choose H^, and consider an in- 
cident electron from the normal side of the junction {x < 0) 
with energy E. For positive excitation energies E > 0, the 
eigenvectors and corresponding momentum of the particles 
read = [1, e'^ 0, 0] V^" p'^ = {E + Ep)/vf, for 

a right-moving electron at angle of incidence 9, while a left- 
moving electron is described by the substitution 9 ^ n — 9. 
If Andreev-reflection takes place, a left-moving hole with 
energy {—E) and angle of reflection 6'a is generated with 
belonging wave-function = [0, 0, 1, e~'^*]^e~'^'''™*'^*'^, 
— [E — Ef:)/vf, where the superscript e (h) denotes 
an electron-like (hole-like) excitation. Since translational in- 
variance in the y-direction holds, the corresponding com- 
ponent of momentum is conserved. This condition allows 
for determination of the Andreev-reflection angle 9a through 



sin 9a = sin 6. One thus infers that there is no Andreev- 
reflection (9a = ±7r/2), and consequently no subgap con- 
ductance, for angles of incidence above the critical angle 

9^ = asin[|£; - Ef\/{E + E^)]. 

On the superconducting side of the system (x > 0), the 
possible wavefunctions for transmission of a right-moving 
quasiparticle with a given excitation energy E > reads 

^-^^ = [u(6l+),u(6l+)e"^^,i;(6l+)e-"^+,t;(6l+)c'(''+-'^+)]T, 
^^iq'cose+x^ ^ ^ _ A^)/vp, and 

^lf^_ = [t;(6i-),f(6'-)e'»^'a(6'-)e-"^^u(0-)e'(«"-'^")]T 
xgiV'cose-^^ ^ _ 1 A^)/vp. The co herence 
factors are given by u{9) = [5(1 + ^/E^~^A{9)^ / E)]^/^, 
v{9) = - ^/E^ - |A(6i)|2/£;)]i/2. Above, we have de- 
fined 6I+ = 91,9- =7r-6lg,ande">'' = e"'A(6'±)/|A(6'±)|. 
The transmission angles for the electron-like 

(ELQ) and hole-like (HLQ) quasiparticles are given by 



p'^ sin 9, i=e,h. Note that for subgap energies 
E < A, there is a small imaginary contribution to the 
wavevector, which leads to exponential damping of the 
wavefunctions inside the superconductor. For clarity, we 
have omitted a common phase factor e'*^"^ which corresponds 
to the conserved momentum in the y-direction. A possible 
Fermi vector mismatch (FVM) between the normal and 
superconducting region is accounted for by allowing for 
E'p 7^ _Ep. The case E'^ ^ E^ corresponds to a heavily 
doped superconducting region, while E'^ = E^ describes 
undoped graphene. Since we are using a mean-field approach 
to describe the superconducting part of the Hamiltonian, it 
is implicitly understood that phase-fluctuations of the order 
parameter must be small. This amounts to imposing the 
restriction fl^ ^/A' > 1, or equivalently, E'^ > A. 

The conductance of the N/I/S junction is given by ifTHl 
G[eV) = Gnf\d9cos9[l-\r{eV,9)\^+P\rA{-eV,9)\% 
where r and ta are the reflection coefficients for normal and 
Andreev reflection, respectively, P = |ph| cos 6'A/(|Pe| cos 9), 
while Gn = /^^/2'16i cos 6l[4cos2 61/(4 cos2 61 + Z"^)] is a 
renormalization constant corresponding to the N/N metallic 
conductance lfl2ll . In this case, we have zero intrinsic bar- 
rier such that Z = 0. We will apply the usual approximation 
|rA(— eV, 9)\ ~ |rA(eV, 9)\, which holds for subgap energies. 
Although it is not valid for energies above the gap, this is of lit- 
tle consequence for the final result, since Andreev reflection is 
suppressed for eV > A. The reflection and transmission co- 
efficients constitute a unitary scattering matrix, a property that 
essentially expresses a conservation of probability. In deriving 
the conductance, we have ensured that the scattering coeffi- 
cients have been normalized by the incoming current through 
the factor P. In order to obtain these coefficients, we make 
use of the boundary conditions '0U=o = i'i\x=0: "^Pilx^d = 
^sU=d, where we have defined the wavefunction in the in- 
sulating region ipi = + t24'- + tai^X + ^4^'-- The 
wavefunctions ip differ from in that the Fermi energy is 
greatly shifted by means of e.g. an external potential, such 
that Ep Ep ~ Vq where Vq is the barrier (equivalent to the 
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role of Z in Ref. The coefficients r and ta may now be 
obtained by using the boundary conditions, but we leave the 
explicit calculations and somewhat cumbersome analytical re- 
sults for a forthcoming paper 

Consider first a N/I/S graphene junction. In the thin-barrier 
limit defined as d — ^ and Vb — * oo with s-wave pairing, Ref. 
Istl reported a 7r-periodicity of the conductance with respect 
to the parameter x — V^dlv-p. In the present study, we do not 
restrict ourselves to isotropic pairing, nor to the thin-barrier 
limit, and show that new physics emerges from the presence 
of a finite-width barrier. We measure the width w of region I in 
units of d/Ap and the potential barrier Vb in units of Ev_. The 
linear dispersion approximation is valid up to ~ 1 eV [80, and 
we will consider typical Fermi energies in graphene of E-p — 
100 meV in the undoped case and a gap A = 1 meV In the 
doped case, we set E'^ = lOEp, and we also fix Vq ~ lOE^ 
in order to operate within the regime of validity of the linear 
dispersion approximation. The undoped situation originally 
refers to the case where the Fermi level is located at the Dirac 
point, although real experimental graphene samples will have 
free carriers, such that Ep is pushed upwards. The doped case 
denotes a large FVM between the N and S region which may 
be induced by chemical doping or by a gate voltage. 
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FIG. 1: Tunneling conductance of N/I/S graphene junction for both 
s-wave and d-wave pairing in the undoped and doped case (see main 
text for parameter values). It is seen that for increasing w, a novel 
oscillatory behaviour of the conductance as a function of voltage is 
present in all cases. 

Consider Fig. [T] where we plot the normalized tunneling 
conductance in the two cases of s-wave and d-wave pairing, 
for both doped and undoped graphene. The most striking new 
feature compared to the thin-barrier limit is the strong oscilla- 
tions in the conductance as a function of eV. We also include 
the thin-barrier limit with % = and x = tt to illustrate the 
TT-periodicity in this limit. For subgap energies, we regain the 
N/S conductance for undoped graphene when x = 0, with 
nearly perfect Andreev reflection. To model the d-wave pair- 
ing, we have used the dj.2_y2 model A{9) — A cos(26' — 2a) 
with a = tt/A. The parameter a effectively models different 
orientations of the gap in k-space with regard to the interface, 
and a = 7r/4 corresponds to perfect formation of ZES in N/S 



metallic junctions. For a = 0, the d-wave spectra are es- 
sentially identical to the s-wave case, since the condition for 
formation of ZES is not fulfilled in this case 101 ■ It is seen that 
in all cases shown in Fig. [T]the conductance exhibits a novel 
oscillatory behavior as a function of applied bias voltage eV 
as the width w of the insulating region becomes much larger 
than the Fermi wavelength, i.e. w ^ Xp. 

The oscillatory behavior of the conductance may be under- 
stood as follows. Non-relativistic free electrons with energy 
E impinging upon a potential barrier Vq are described by an 
expontentially decreasing non-oscillatory wavefunction c'''^ 
inside the barrier region if E < Vq, since the dispersion es- 
sentially is fc ~ \/E — Vq. Relativistic free electrons, on the 
other hand, have a dispersion A; ~ {E — Vq), such that the 
corresponding wavefunctions do not decay inside the barrier 
region. Instead, the transmittance of the junction will display 
an oscillatory behavior as a function of the energy of inci- 
dence E. In general, a kinetic energy given by ~ k" will 
lead to a complex momentum k ^ {E — Vq)^^" inside the 
tunneling region, and hence damped oscillatory behavior of 
the wave function. Relativistic massless fermions are unique 
in the sense that only in this case {a = 1) is the momentum 
purely real. Hence, the undamped oscillatory behavior at sub- 
gap energies appears as a direct manifestation of the relativis- 
tic low-energy Dirac fermions in the problem. This observa- 
tion is also linked to the so-called Klein paradox which occurs 
for electrons with such a relativistic dispersion relation, which 
has been theoretically studied in normal graphene [2]. 

We next discuss why the illustrated conductance spectra 
are different for s-wave and d-wave symmetry, in addition 
to comparing the doped and undoped case. The doping level 
may be considered as an effective FVM, acting as a source 
of normal reflection in the scattering processes. This is why 
the subgap conductance at thin barrier limit is reduced in the 
doped case. Moving away from the thin barrier limit, it is 
seen that oscillations emerge in the conductance spectra. For 
s-wave pairing, the amplitude of the oscillations is larger in 
the doped case than in the undoped case, and the period of os- 
cillations remains the same. This period depends on w, while 
the amplitude of the oscillations is governed by the wavevec- 
tors in the regions I and S. The maximum value of the oscil- 
lations occurs when 2w equals an integer number of wave- 
lengths, corresponding to a constructive interference between 
the scattered waves. Physically, the amplitude-dependence of 
the oscillations on doping originates with the fact that doping 
effectively acts as an increase in barrier strength. By mak- 
ing Vq larger, one introduces a stronger source of normal re- 
flection. When the resonance condition for the oscillations is 
not met, the barrier reflects the incoming particles more effi- 
ciently. This is also the reason why increasing Vq directly and 
increasing Ep has the same effect on the spectra. 

We now turn to the difference between the s-wave and d- 
wave for the undoped case. It is seen that the conductance is 
reduced in the d-wave case compared to the s-wave case. One 
may understand the reduction in subgap conductance in the 
undoped case as a consequence of tunneling into the nodes of 
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the gap, which is not present in the s-wave case. Hence, An- 
dreev reflection which significantly contributes to the conduc- 
tance, is reduced in the d-wave case compared to the s-wave 
case. Moreover, we see that a ZBCP is formed in the doped 
case, equivalent to a stronger barrier, and this is interpreted as 
the usual formation of ZES leading to a transmission at zero 
bias with a sharp drop for increasing voltage. 
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FIG. 2: Conductance spectra for a doped N/S graphene junction in 
the d-wave case, varying the orientation angle a. From right to left, 
the curves vary in the interval a £ [0,7r/4] in steps of 7r/20. The 
inset shows the undoped case for a = 7r/4. 



a = 7r/4, where the deviation from perfect Andreev reflection 
for eV < A is due to tunneling into the nodes of the gap. 

In summary, we have studied coherent quantum transport in 
N/S and N/I/S graphene junctions, investigating also the role 
of d-wave pairing symmetry on the tunneling conductance. 
We report a new oscillatory behaviour of the conductance as 
a function of bias voltage for insulating regions that satisfy 
d ^ \, which is present both for s- and d-wave pairing. In 
the latter case, we have also studied the conductance of an N/S 
junction and find very distinct behaviour from metallic N/S 
junctions: a rotation of a is accompanied by a progressive 
shift of the peak in the conductance. All of our predictions 
should be easily experimentally observable. 
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Finally, we briefly investigate how the conductance spec- 
tra of a N/S graphene junction (without the insulating region) 
change when going from a s-wave to a d-wave order parame- 
ter in the superconducting part of the system. Consider Fig. |2] 
for the case of doped graphene, where we plot the conductance 
to see how it evolves upon a rotation of the gap. The behaviour 
is quite distinct from that encountered in a N/S metallic junc- 
tion. From Fig. |2] we see that the peak of the conductance 
shifts from eV = A to progressively lower values as a in- 
creases from to 7r/4, but only for a very close to 7r/4 is 
a ZBCP present. This is different from what is observed in 
metallic N/S junctions, where the formation of a ZBCP starts 
immediately as one moves away from a = in the presence 
of a FVM, corresponding to the doped case here. In Fig. |2]the 
conductance spectra actually mimicks a lower value of the gap 
than what is the case, if one were to infer the gap magnitude 
from the position of the singularity in the spectra. This should 
be an easily observable feature in experiments, and provides 
a direct way of testing our theory. For undoped graphene, we 
found very little difference in the conductance spectra upon 
varying a. The inset of Fig. |2]illustrates the undoped case for 
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